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ATTACHMENT  II 


ABSTRACT 


'  .„^The  parameters  in  a  linear  regression  model  can  be  estimated  by 
minimizing  the  sum  of  the  absolute  residuals  (L^  estimation)  Instead 
of  the  more  classical  approach  of  minimizing  the  sum  of  squared 
residuals  (least  squares  estimation).  In  addition  to  other  nice 
properties  estimators  are  less  sensitive  to  outliers  than  least 
squares  estimators.  This  paper  describes  a  linear  programming 
algorithm  and  computer  program  for  obtaining  estimators  and 
estimates  of  their  covariances  when  the  regression  parameters  are 
restricted  to  satisfy  specified  linear  constraints.  These  estimated 
covariances  are  the  new  feature  in  this  work  and  are  an  extremely 
important  ingredient  in  hypothesis  tests  and  confidence  interval 
construction.^  Technical  Report  64  describes  a  similar  procedure  for 
obtaining  unbiased  estimators  when  there  are  no  constraints  on 
the  parameters. 
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Restricted  Estimators  and  their  Covariances 

1.  An  Introduction  to  MR  ■  A 

Consider  the  linear  regression  model  in  the  form 

y  -  X0  +  e  (1) 

where  y  is  a  vector  of  n  observations,  X  is  an  n  x  p  matrix  of  rank  p  of 

known  constants,  g  is  a  vector  of  p  unknown  parameters  and  c  is  a  vector 

of  independent  random  variables  (noise)  symmetrically  distributed  with 

2 

mean  zero  and  variance  a  .  The  estimation  of  g  can  be  obtained  under 
several  different  optimality  criteria.  For  g  unrestricted,  the  classical 
least  squares  estimator,  g,  where 

g  =  (XTX)_1  XTY  ,  (2) 

has  the  smallest  variance  among  the  class  of  unbiased  linear  functions  of 
y.  However,  the  least  squares  estimator  is  extremely  sensitive  to  large 
values  of  |e|,  outliers,  particularly  when  the  sample  size,  n,  is  small 
•relative  to  p,  say  n  <_  2(p+l).  In  addition,  the  least  squares  estimator 
does  not  have  the  flexibility  of  allowing  restrictions  to  be  placed  on  g. 
These  two  drawbacks  suggest  that  an  optimality  criteria  other  than  least 
squares  be  considered.  Several  authors  (Barrodale  (1968),  Charnes  and 
Cooper  (1964),  Gentle,  Kennedy  and  Sposlto  (1977),  Harris  (1950), 

Harter  (1974),  Rice  and  White  (1964),  Taylor  (1973)  have  suggested  that 

n 

z  !  y4  -  x  g|  (3) 

i-i  1  1 

should  be  minimized  with  respect  to  g  where  y^  is  the  i-th  observation  and 

A 

is  the  i-th  row  of  X.  The  estimator,  g,  which  minimizes  the  sum  of  the 
absolute  residuals  is  often  called  the  estimator. 

The  estimate  is  not  necessarily  unique  since  it  depends  upon  the 
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method  of  computation.  However,  conventional  linear  programming  algorithms 
can  calculate  the  estimate  with  m  added  linear  restrictions  on  0  in  the 
forms, 

A0  <_  b  or  (4) 

Ag  ^  b  or 

A0  =  b 

where  A  is  an  m  x  p  matrix  of  known  coefficients  for  g  and  b  is  a  vector 
of  m  constants. 

The  computer  program  MR.  A  implements  an  algorithm  for  Minimizing  the 
sum  of  the  Absolute  Residuals  where  0  is  restricted  by  any  of  the  relation- 
ships  given  in  (4).  The  algorithm  uses  the  least  squares  estimator  0  as  an 
initial  unbiased  antisymmetrical  estimator  of  0  where 

0  -  0  (e)  =  -  [0  -  0  (-£)]  (5) 

in  order  to  provide  efficiency  in  computation  as  suggested  by  Barrodale 
and  Roberts  (1973). 

The  estimation  of  the  unrestricted  estimator  is  only  a  first  step 
in  the  completion  of  the  task  of  estimating  0  in  model  (1)  subject  to  re¬ 
strictions  in  (4).  The  second  step  is  the  estimation  of  the  covariances 

A 

of  0.  This  estimation  would  usually  be  required  for  any  further  develop¬ 
ments  in  hypothesis  testing  or  confidence  intervals.  MR.  A  contains  a 
mini-Monte  Carlo  procedure  similar  to  that  used  in  the  algorithm  MRS.  A 
(Book,  Booker,  Hartley,  Sielken  1980)  for  estimating  the  covariance  of 
the  restricted  estimator.  This  feature  sets  MR.  A  apart  from  other 
restricted  procedures. 

2.  Computational  Procedure 

The  problem  of  minimizing  the  sum  of  the  absolute  residuals  can  be 
formulated  as  follows: 

n 

■ln  i!x  <ru  +  'a1 


(6) 


subject  to 


yi  "  Xi0  "  (rli  "  r2i)*  1  "  1 . . 


rli»  r2i  -  °» 


(7) 


<  > 


V 


where  the  i~th  residual,  r^,  is  expressed  as  the  difference  between  two 
positive  components,  r^  and  Prom°te  the  unbiasedness  of  the 

resulting  estimator  and  to  increase  the 

efficiency  of  computation  in  the  algorithm,  the  problem  is  reformulated 

following  Hartley  and  Sielken  (1973)  by  introducing  the  antisymmetrical 

* 

least  squares  estimator,  3q  =6,  and  transforming  (7)  to 


y±  “  -  X±  (3  -  8q)  a  ~  ^2i*  ^  •••»  ^ 


(8) 


Ai(e  -  B0)  > 


>  bi  -  Vo 


Then,  using 


8  -  @(1)  -  e(2)  , 

•o  ■  -  42) 

with  3(1),  8^2),  Bq2)  >  0  in  (8)  yields 

*i (B<1)  +'bo2>)  -  xi  (b<2)  +  BoU)  +  rn  -  r2i  -  n  -  v 

(S(1)  +  8j2))  -  Aj  (b(2)  +  ’  bt  -  A160 

or  equivalently 


0 

(9) 


XiBl  "  Xi®2  +  rli  “  r2i  “  yi  ~  XiB0 


AiBl  “  AiB2 


{•}  “•  ‘ v 


for  i  **  1,  n  where 


B1  -  8(1)  +  >_  0 


B2  =  8(2)  +  8qX)  >  0  . 


Now,  in  order  to  promote  unbiasedness  of  8  and  to  partially  compensate  for 
any  idiosyncrasies  in  the  particular  linear  programming  algorithm  used  to 
solve  the  problem  in  (6),  (7),  (9),  (10),  the  problem  is  considered  in  two 
equivalent  forms  and  P2  with  MR.  A  randomly  selecting  either  P^  or  P2 
with  probability  1/2.  In  matrix  notation,  the  problems  P^  and  P2  are 


as  follows: 


V  ”ln  (rll  +  r2i> 


subject  to 


X  -X 


A  -A 


1  -1 


0  0 


r  i 

i  ei  i 


“2  * 


r  1 
I  y  -  XBo  I 


b  -  A8. 


Bj  ,  B2  ,  r^  ,  r2  _>  0; 


P2:  min  ^  <ru  + 


subject  to 
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B 


1 


* 


>  0. 


After  HR.  A  selects  and  solves  either  or  P the  estimator  is  B 
where 

6  -  Bi  -  b2  +  BQ.  (11) 

2 

In  cases  where  the  sample  size  is  large,  say  n  >>  (p+1)  ,  the  sample 

A 

could  be  randomly  divided  into  G  groups,  then  B  estimated  for  each  group 

S 

A 

g  separately,  and  the  covariance  of  6  estimated  from  the  sample  covariance 

of  the  0  's. 

g 

MR.  A  estimates  the  covariance  of  B  for  any  sample  size  using  a  mini- 
Monte  Carlo  procedure.  A  Monte-Carlo  estimate  is  obtained  by  essentially 
generating  several  sets  of  n  y's,  finding  the  estimate  for  each  set,  and 
computing  the  sample  covariance.  There  are  two  difficulties  with  generating 
the  y’s;  namely. 


(1)  B  is  unknown,  and 
2 

•  (2)  o  is  unknown. 

To  overcome  the  first  difficulty,  the  sum  of  the  absolute  residuals  is  ex¬ 
pressed  as 


n 

Z 

i-1 


n 

Z 

i-1 


V 


+  Ej 


-v 


-  Z  I  €  -  X.(B-fS)! 
i-1  1 


n 


(12) 


-  Z  |  e  -  X  6B  | 
i-1  1 

A 

using  (1)  and  66  -  B-fJ.  Note  that  the  covariance  of  the  estimator  of 

A 

66  is  the  same  as  the  covariance  of  B»  Thus,  only  sets  of  n  c's  need  be 

i 

generated.  To  deal  with  the  second  difficulty  (the  unknown  o  ),  note  that 
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Z  I  £  -  X  66 

i-1 


“  JL  ^ 

I  Z  I  e.  -  X  68 
i=l 


and  the  restrictions  on  8  in  (A)  become 


A8  ~  A8 


A(68) 


G) 

(!) 


b  -  AB 


b  -  A6 


A(68  )  <  >  >  b  -  A6 


where  66  =  (8-8) /o  and  where  =  e^/c  and  is  symmetrically  distributed 

with  mean  0  and  variance  1.  Furthermore, 

*  *  2  *  * 

Covariance  (8)  =  Covariance  (68)  *  o  Covariance  (68  ).  (15) 

* 

The  mini-Monte  Carlo  procedure  generates  K  sets  of  n  e  's;  finds 

*  *  * 
the  K  estimates  68  ;  and  estimates  the  covariance  of  8  by  o  times  the 

sample  covariance  of  68  where  a  is  calculated  from 


o  -  I  |  y  -  X  6  | / (n-p) 
i-1  x 

- - 

k  jj  jj  i  4  -  vK  i  '<-»> 


This  estimator  of  a  reflects  both  the  variability  in  y^  -  X^B  and  the 
variability  associated  with  the  linear  programming  algorithm. 

In  the  estimation  of  each  68^,  an  additional  problem  arises  from 
the  right  hand  side  of  (14)  being  a  function  of  8  and  o.  This  problem  is 
overcome  by  estimating  8  using  (11)  and  estimating  a  using  (16)  for 
each  sample,  k  >  2,  as  it  is  generated.  For  the  first  sample  (k  -  1), 
o  is  estimated  by  the  numerator  of  (16).  "his  numerator  has  been  shown 
to  provide  a  reasonable  alternative  ast  tator  of  o  [Book,  Booker,  Hartley, 


-7- 


MR.  A  allows  the  user  to  generate  the  e  's  from  the  uniform, 
normal  or  double  exponential  distributions.  These  distributions  were 
selected  in  order  to  represent  short,  medium  and  long  tailed  distributions, 
respectively.  These  three  distributions  are  also  interesting  because 
maximum  likelihood  estimation  corresponds  to  minimizing  the  maximum  absolute 
residual  if  the  e's  are  uniform,  minimizing  the  sum  of  squared  residuals  if 
the  e's  are  normal,  and  minimizing  the  sum  of  the  absolute  residuals  if  the 
e's  are  double  exponential. 

MR.  A  also  allows  the  user  to  assign  weight  coefficients  to  the 
residuals,  changing  (6)  to 

n 

min  Z  w  (r  +  r,  ) 
i=l  1  11 

where  is  the  weigh!  coefficient  for  the  i-lh  lesiuuel,  r^.  If  the 
w^'s  are  not  all  equal  to  one,  then  the  objective  function  in  the  mini- 
Monte  Carlo  study  becomes 

^  wi  I  ei  -  V**  i 

i=l  1  1 

and  the  estimate  of  a  also  reflects  the  w^,  namely, 

»  n 

o  -  Z  w  |  y  -  X  6  |  /(n-p) 
i-1  1  1 

k  X  “i  i  4  -  vs»\  i  '<-»> 

k*l  i=l 

3.  MR.  A  Users’  Guide  and  Sample  Problem 

MR.  A  consists  of  a  main  program  and  subroutines.  The  functions  of 
these  components  are  given  in  Table  1  at  the  end  of  this  section. 

The  input  to  the  computer  program  consists  of  five  basic  card  groups. 

Input  Instructions  are  also  briefly  documented  in  the  program  comment  statements. 
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The  first  card  or  card  image  contains  the  ancillary  statistics  re¬ 
quired  for  the  specific  problem  as  follows: 


:  Card: 

Card  Column 

Variable  Name 

Description 

1-3 

IP 

=  p  =  number  of  beta  parameters  (format 

13;  i.e.,  a  3  digit  integer,  right  justified) 

4-6 

NOBS 

=  n  »  number  of  observations  (format  13) 

7-9 

NC 

=  m  =  number  of  constraints  (format  13) 

11-20 

NSEED 

Ten  digit  random  number  <  2147483647 

(format  HO) 

21-25 

ISAM 

=  K  =  number  of  samples  for  the  mini- 

Monte  Carlo  study  (format  15) 

27 

IuRIT0 

a  1  if  ^robl^ni  irnut  is  to  b n  nrintc - 

*  0  otherwise 

29 

IWRIT1 

=  1  if  the  main  results  of  the 

estimation  are  to  be  printed 

*  0  otherwise 

31 

IWRIT2 

=  1  if  the  intermediate  results  of  the 

estimation  are  to  be  printed 

**  0  otherwise 

33 

IWRIT3 

*  1  if  the  basis  inverse  is  to  be  printed 

*  0  otherwise 

35 

1 WRIT 4 

*  1  if  the  main  results  of  the  mini- 

Monte  Carlo  study  are  to  be  printed 

■  0  otherwise 

37 

IWRIT5 

■  1  if  the  intermediate  results  of  the 

mini-Monte  Carlo  study  are  to  be  printed 

0  otherwise 
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L 


Card  Column 


Variable  Name 


IWRIT6 


Description 


1  If  the  Intermediate  results  In  the 


determination  of  the  covariance  of  B 


are  to  be  printed 


0  otherwise 


10PTN 


1  11  the  e  's  are  to  be  normally  dis¬ 


tributed 


2  if  the  e  's  are  to  be  double  ex¬ 


ponentially  distributed 

* 

*  3  If  the  e  's  are  to  be  uniformly 


distributed. 


1  If  the  residuals  are  assigned  weights. 


0  if  the  residuals  are  not  weighted. 


The  remaining  card  input  instructions  are  as  follows: 


Card  Type 


Variable  Name 


Description 


Second  card  group:  ITYPE^  i*=l,  ....  NC  The  type  of  the  i-th  constraint  on  B 

■  0  if  the  constraint  is  <_  type 

■  1  if  the  constraint  is  >_  type 

*  2  if  the  constraint  is  =  type 

(Format  12) 

Third  card  group:  y^,  i»l,  NOBS  The  observations  [format  (10F8.3); 

i.e.,  ten  eight  digit  numbers  with 
either  a  decimal  point  included  or 
last  five  digits  are  assumed  to  be 
to  the  right  of  a  supplied  decimal 
point] 


Variable  Name 


Description 


Card  Type  Variable  Name 

x^,  i-1,  NOBS  The  matrix  of  beta  coefficients. 


J-l.  ....  IP 


Fourth  card  group:  A^,  i*l,  NC 


J-l ,  •  •  •  ,  IP 


bt,  1-1, 


,,  NC 


read  in  by  rows  [format  (10F8.3)]. 
The  matrix  of  beta  coefficients 
in  the  constraints  on  beta,  read  in 
by  rows  [format  (10F8.3)). 

The  right-hand  sides  of  the  beta 
constraints  [format  (10F8.3)]. 


Fifth  card  group 


w^,  i=l,  NOBS  The  weights  assigned  to  the  re¬ 


siduals  [format  (10F8.3)]. 
(Included  only  if  IWT  =  1)  . 

The  user  then  supplies  a  title  card  of  length  80  spaces  or  less  fol¬ 
lowing  the  fifth  card  group. 

The  size  of  the  problem  which  MR.  A  can  solve  is  limited  only  by  the 
current  dimensioit  sLatciiicntS.  These  restrict  the  number  of  observations  to 
20  (n  20),  the  number  of  parameters  to  10  (p  <_  10),  the  number  of  con¬ 
straints  to  10  (m  10) ,  and  the  number  of  samples  in  the  mini-Monte  Carlo 
study  to  100  (K  <_  100) .  However,  expansion  can  easily  be  accomplished  by 
increasing  these  dimensions  in  the  dimension  statements  as  documented  in 
the  program. 

MR.  A  is  written  in  Fortran  IV  language  and  is  compatible  with 
Fortran  G  and  H  and  WATFIV  language  compilers.  The  program  uses  double 
precision  arithmatic. 

MR.  A  has  been  tested  on  several  problems  on  an  AMDAHL  470  V6  and 
should  be  compatible  with  all  IBM  compilers.  MR.  A  will  execute  small 
problems  such  as  n  -  5,  p  -  3,  K  -  6,  m  «*  3  in  less  than  two  seconds. 
Problems  of  the  size  n  -  20,  p  -  2,  K  “  30,  m  -  3  take  up  to  a  minute  of 
execution  time.  A  sample  problem  with  sample  Input  and  output  is  given 
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in  Appendix  A  and  B  respectively.  The  program  listing  of  MR.  A  is  given 
in  Appendix  C. 


'i 
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Table  1 


Component 

MAIN 


DUAL 

MCRHS 

SELECT 

RAND 

RHANDS 

BTILDE 

INVERT 

NORMAL 


DOUBLE 


UNIFRM 


Components  of  MR.  A  and  their  Functions 

Func  t Ion 

Reads  data  and  generates  output. 

A 

Performs  L^  estimation  of  B. 

•*»  A 

Determines  a  and  the  estimated  covariance  of  B. 

Carries  out  the  mini -Monte  Carlo  study. 

Forms  the  right  hand  sides  of  (14)  for  the  mini- 
Monte  Carlo  study. 

Selects  P^  or  P^  for  solution. 

Generates  random  uniform  variables  ranging  from  0  to  1. 
Forms  the  right-hand  sides  of  (9). 

Constructs  the  least  squares  estimate  of  B;  i.e.,  Bq. 
Inverts  an  n  x  n  matrix. 

Generates  a  vector  of  normally  distributed  random 
variables  with  mean  0  and  variance  1  for  the  mlni- 
Monte  Carlo  study. 

Generates  a  vector  of  double  exponentially  distributed 
random  variables  with  mean  0  and  variance  1  for  the 
mini-Monte  Carlo  study. 

Generates  a  vector  of  uniformly  distributed  random 
variables  with  mean  0  and  variance  1  for  the  mini- 


Monte  Carlo  study. 
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APPENDIX  A.  SAMPLE  INPUTS 


253  1872539680  20  0  O  0  0 

2  10 

4.257  -1.983  0.024  -3.180 

1.000  1.120 

1.000  1.950 

1.000  3.020 

1.000  5.430 

1.000  6.590 

1.000  3.000  -1.000 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MR. 


0  0  0  1  0 
-3.986 


A.  WITH  NO  OPTIONAL  PRINTOUTS. 


253  1872539680  20  1  1  1  1 

2  10 

4.257  -1.983  0.024  -3.180 

1.000  1.120 
1.000  1.950 

1 .000  ‘  3.020 

1.000  5.430 

1.000  6 .590 

1.000  3.000  -1.000 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MR, 


1  I  1  1  0 

-3.986 


A.  ttITH  ALL  OPTIONAL  PRINTOUTS, 
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APPENDIX  B. 


SAMPLE  OUTPUTS 
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MR.  A.  1 

HINIHIZES  THE  ABSOLUTE  RESIDUALS. 

THIS  PROGRAM  ESTIMATES  A  LINEAR  REGRESSION  BY  MINIMIZING  (Sample  output  with  no  optional  printouts) 
THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  -  LI  ESTIMATION. 

IN  ADDITION.  A  MINI-MONTE  CARLO  SIMULATION  GENERATES  AN 
ESTIMATED  COVARIANCE  MATRIX  FOR  THE  ESTIMATED 
REGRESSION  PARAMETERS. 

IN  ADDITION.  LINEAR  CONSTRAINTS  MAT 
BE  PLACED  ON  THE  PARAMETER  VECTOR.  BETA, 

THESE  CONSTRAINTS  MAY  BE  EQUALITY  OR  INEQUALITY  CONTSTRAINTS. 

THE  FOLLOWING  PROCEDURE  DEVELOPED  BY  J 
D.N.  BOOK 
J.B.  BOOKER 
H.O.  HARTLEY 
R.L.  SIELKEN.  JR. 

INSTITUTE  OF  STATISTICS 
TEXAS  A  1  M  UNIVERSTIY 
COLLEGE  STATION.  TEXAS  77843 

INQUIRIES  AND  COMMENTS  SHOULD  BE  ADDRESSED  TO*. 

ROBERT  L.  SIELKEN.  JR. 

THE  SUPPORT  OF  THE  OFFICE  OF  NAVAL  RESEARCH  IS  GRATEFULLY  ACKNOWLEDGED. 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MR.  A.  WITH  NO  OPTIONAL  PRINTOUTS. 

NUMBER  OF  OBSERVATIONS  =  5 

NUMBER  OF  PARAMETERS  =  2 

USER  SUPPLIED  RANDOM  INTEGER,  NSEED  =  1872539680 

SAMPLE  SIZE  FOR  THE  MINI  MONTE  CARLO  STUDY  =  20 

NUMBER  OF  BETA  CONSTRAINTS  =  3 

NUMBER  OF  <  OR  =  BETA  CONSTRAINTS  =  1 

NUMBER  OF  >  OR  =  BETA  CONSTRAINTS  =  1 

NUMBER  OF  =  BETA  CONSTRAINTS  =  1 

THE  AUXILIARY  LEAST  SQUARES  ESTIMATE,  BETAO,  OF  THE  REGRESSION  PARAMETER  VECTOR,  BETA 
LEAST  SQUARES  ESTIMATE  OF  BETA(  1)  =  3.227002 

LEAST  SQUARES  ESTIMATE  OF  BETA!  2)  *  -1.159747 

MR.A'S  ANSWER  i  THE  ESTIMATE  OF  THE  REGRESSION  PARAMETER  VECTOR  WHICH  MINIMIZES  THE  SUM  OF  THE  ABSOLUTE  RESIDUALS 


LI  ESTIMATE  OF  BETA!  1)  =  3.000000 

LI  ESTIMATE  OF  BETA!  2)  =  -2.000000 

THE  RESIDUALS  1 


THE 

1-TH  RESIDUAL  = 

3.497000 

THE 

2-TH  RESIDUAL  = 

-1.083000 

THE 

3-TH  RESIDUAL  = 

3.064000 

THE 

4-TH  RESIDUAL  = 

4.680000 

THE 

5-TH  RESIDUAL  = 

6.194000 

THE  SUM  OF  THE  ABSOLUTE  RESIDUALS 

-  18.518000 

THE  MAXIMUM  ABSOLUTE  RESIDUAL  = 

6.194000 

MAIN  RESULTS  OF  THE  MINI  MONTE  CARLO  STUDY 


ESTIMATED  VALUE  OF  SIGMA 


ESTIMATED 


3.56543? 

-0.970588 


OF  Tl 


2*817978 

PARAMETER  VECTOR,  BETA 

J.970588 
0.435623 


MINIMIZES  THE  ABSOLUTE  RESIDUALS 


THIS  PROGRAM  ESTIMATES  A  LINEAR  REGRESSION  BY  MINIMIZING 
THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  -  LI  ESTIMATION.  (Sample  out 

IN  ADDITION.  A  MINI-MONTE  CARLO  SIMULATION  GENERATES  AN 
ESTIMATED  COVARIANCE  MATRIX  FOR  THE  ESTIMATED 
REGRESSION  PARAMETERS. 

IN  ADDITION.  LINEAR  CONSTRAINTS  MAY 
BE  PLACED  ON  THE  PARAMETER  VECTOR.  BETA. 

THESE  CONSTRAINTS  MAY  BE  EQUALITY  OR  INEQUALITY  CONTSTRAINTS. 

THE  FOLLOWING  PROCEDURE  DEVELOPED  BY  S 
D.N.  BOOK 
J.B.  BOOKER 
H.O.  HARTLEY 
R.L.  SIELKEN.  JR. 

INSTITUTE  OF  STATISTICS 
TEXAS  A  t  M  UNI VERST I Y 
COLLEGE  STATION.  TEXAS  77843 

INQUIRIES  AND  COMMENTS  SHOULD  BE  ADDRESSED  TOi 
ROBERT  L.  SIELKEN.  JR. 

THE  SUPPORT  OF  THE  OFFICE  OF  NAVAL  RESEARCH  IS  GRATEFULLY  ACKNOWLEDGED. 

THIS  IS  AN  EXAMPLE  PROBLEM  OF  MR.  A.  WITH  ALL  OPTIONAL  PRINTOUTS. 

NUMBER  OF  OBSERVATIONS  =  5 

NUMBER  OF  PARAMETERS  =  2 

USER  SUPPLIED  RANDOH  INTEGER,  NSEED  =  1872539680 

SAMPLE  SIZE  FOR  THE  MINI  MONTE  CARLO  STUDY  =  20 

NUMBER  OF  BETA  CONSTRAINTS  =  3 

NUMBER  OF  <  OR  =  BETA  CONSTRAINTS  =  1 

NUMBER  OF  >  OR  =  BETA  CONSTRAINTS  =  1 

NUMBER  OF  =  BETA  CONSTRAINTS  =  1 

THE  LINEAR  REGRESSION  IS  Y  =  XlBETA  I  EPSILON 
WHERE 

Y  CONTAINS  THE  OBSERVATIONS, 

BETA  IS  A  VECTOR  CONTAINING  THE  REGRESSION  PARAMETERS, 

THE  I-TH  ROW  OF  X  CONTAINS  THE  COEFFICIENTS  OF  BETA 

CORRESPONDING  TO  THE  I-TH  OBSERVATION,  AND 

EPSILON  IS  A  RANDOM  VARIABLE  WITH  MEAN  ZERO  AND  VARIANCE 

SIGMA-SQUARED  REPRESENTING  RANDOH  VARIABILITY  FROM 

XtBETA. 

THE  Y  VECTOR 

4.25700 

-1.98300 

0.02400 

-3.18000 

-3.98600 

THE  X  MATRIX 

1.0  1.1 
1.0  2.0 
1.0  3.0 

1.0  5.4 

1.0  6*6 

THE  CONSTRAINTS  PLACED  ON  BETA  OF  THE  FORM  A4BETA  =  B 
WHERE  THE  CONSTRAINTS  MAY  BE  EQUALITY  OR  INEQUALITY  CONSTRAINTS. 

A  =  THE  HATRIX  OF  COEFFICIENTS  WHOSE  I-TH  ROW  CORRESPONDS  TO  THE  I-TH  INEQUALITY. 
THE  RIGHT  HAND  SIDES  OF  THE  BETA  CONSTRAINTS,  THE  VECTOR  B 
1.00000 
3.00000 
-1.00000 

THE  COEFFICIENTS  FOR  THE  BETA  CONSTRAINTS,  THE  A  MATRIX 
1.0000  1.0000 
1.0000  0.0000 
0.0000  1.0000 
THE  INVERSE  OF  XTX 

0,810298  -0.168497 

-0.168497  0.046521 

THE  VALUE  OF  BETAO  TO  COMPUTE  RHS 
3.227002  -1.159747 

THE  RIGHT  HAND  SIDES,  Y-XIBETAO 

2.32891  -2.94850  0.29943  -0.10958  0.42973 

THE  RIGHT  HAND  SIDES,  B-AIBETAO 
-1.06726  -0.22700  0.15975 

THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  LINEAR  PROGRAMMING  PROBLEM  AS  IT  WAS  CREATED 
THE  OBJECTIVE  FUNCTION  COEFFICIENTS 


(Sample  output  with  optional  printouts) 


ffl 


THE  RIGHT-HAND  SIDES 

BRHS< 

1)  = 

2.32891 

PRHS( 

2)  = 

-2.94850 

BRHS( 

3)  = 

0.29943 

BRHSt 

A)  = 

-0.10958 

BRHSt 

5)  = 

0.42973 

FRHS( 

6)  = 

-1.06726 

BRHSt 

7)  = 

-0.22700 

BRHSt 

8)  = 

0.15975 

1.00000 

1.00000 

-0.00000 


THE  CONSTRAINT  MATRIX  A 

-0.00000 

0  1.00000 

0  1.00000 


0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

1.00000 

-1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 


0.00000 

0.00000 

0.00000 

-1.00000 

0.00000 

0.00000 

-1.00000 

0.00000 

0.00000 

-1.00000 

1.00000 

-1.00000 

-1.00000 

0.00000 

0.00000 

-1.00000 

0,00000 

0.00000 

-1.00000 

0.00000 

0.00000 

-1,00000 

0.00000 

0.00000 


-0.00000 

0.00000  0.00000  0.00000  0,00000  0.00000 

0.00000  0.00000  0.00000  0.00000  0.00000 

1.00000 

THE  INITIAL  BASIS  INVERSE 

Tiff  1  -TH  ROUS  0.10000000D  01  0,000000000  00  -0.10000000B 

0.10000000D  01  -0.10000000D  01  0,000000000  00  O.OOOOOOOOD  00 

THE  2  -TH  ROUS  O.OOOOOOOOD  00  0.10000000D  01  O.OOOOOOOOD 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  0.10000000D  01  O.OOOOOOOOD  00 

THE  3  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  0.10000000D 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

THE  4  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

THE  5  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

THE  6  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

0.10000000D  01  0.000000000  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

TIC  7  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 


-0.00000 

1.00000 

1.00000 

-0.00000 

1.00000 

1.00000 

-0.00000 

1.00000 

-0.00000 

1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

-1.12000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

1.12000 

-1.00000 

0.00000 

1.00000 

-1.95000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

1.95000 

0.00000 

0,00000 

0.00000 

-3.02000 

0.00000 

0.00000 

1.00000 

0.00000 

0.00000 

3.02000 

0.00000 

0.00000 

0.00000 

-5.43000 

1.00000 

-1.00000 

1,00000 

0.00000 

0.00000 

5.43000 

0.00000 

0.00000 

0.00000 

-6.59000 

0.00000 

0.00000 

1.00000 

1.00000 

-1.00000 

6.59000 

0.00000 

0,00000 

0.00000 

-1.00000 

0.00000 

0.00000 

1,00000 

0.00000 

0,00000 

1.00000 

0.00000 

0.00000 

0.00000 

-0.00000 

0.00000 

0.00000 

1.00000 

0,00000 

0.00000 

0.00000 

0.00000 

-1.00000 

0.00000 

-1.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0,00000 

1.00000 

0.00000 

0.00000 

0.00000 

O.OOOOOOOOD  00  0. 10000000D  01  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

THE  8  -TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 
O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0.10000000D  01  O.OOOOOOOOD  00 


01  0.10000000B  01 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
01  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -0.10000000D  01 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 


-0.10000000D  01 
O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
0.10000000D  01 
O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
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THE  9  -TH  R0«:  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  t.OOOOOOOOP  00 

0 >000000000  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0. 100000000  01  O.OOOOOOOOD  00 

THE 10  -TH  ROW:  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  0.10000000D  01 

THE  VARIABLES  INITIALLY  IN  THE  AUGMENTED  BASIS  MATRIX 
0.50.  5.11,  7,13,  9,51,16,15, 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 
XBt  1)  =  -0.61161D  01 

XB(  2)  =  -0.10673D  01 

XB<  3)  =  0.23289D  01 

XB(  4)  =  0.29485D  01 

XB(  5)  =  0.29943H  00 

XB(  6)  =  0.10958D  00 

XB(  7)  =  0.42973D  00 

XBt  8)  =  0.10673D  01 

XB<  9)  =  0.22700D  00 

XBt  10)  =  0.15975D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.1067255616B  01 

PHASE  1  NOW  BEGINS. 

PHASE  1  1  ITERATION  1 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  0.1000000000D  01 

THE  2-TH  REDUCED  COST  =  0.1000000000D  01 

THE  3-TH  REDUCED  COST  =  -0.1000000000D  01 

THE  4-TH  REDUCED  COST  =  -0.10000000008  01 

THE  5-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  7-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  10-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  1 1-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  12-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  13-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  H-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  15-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

X<  1)  IS  THE  NON-BASIC  VARIABLE  WITH  THE  SMALLEST  INDEX 
WHICH  HAS  THE  LARGEST  POSITIVE  REDUCED  COST. 

THEREFORE  X<  1)  WILL  ENTER  THE  BASIS. 

THE  COMPONENTS  OF  Y<-,  1) 

Y(  3,  1)  =  -0.1000000000D  01 

Y(  4,  1)  =  0.1000000000D  01 

Y<  5,  1)  =  -0.1000000000D  01 

Y<  6,  1)  =  0.1000000000D  01 

Y<  7,  1)  =  -0.1000000000D  01 

Y(  8,  1)  =  0.1000000000D  01 

Y(  9,  1)  =  0.1000000000D  01 

Y<  10,  1)  =  O.OOOOOOOOOOD  00 

THE  VECTOR  TO  LEAVE  THE  BASIS  IS  NOW  DETERMINED. 

THE  XB( I )/Y ( I »  1)  RATIOS  WITH  Y(I,  1)  >  0.  ARE! 

0.29485D  01=  XBt  4)/Y(  4,  1) 

0.10958D  00=  XB(  6)/Y(  6,  1) 

0.10673D  01=  XB(  8)/Y(  8,  1) 

0.22700D  00=  XB<  9)/Y(  9,  1) 

THE  MINIMUM  RATIO  IS  XBt  6)/Y<  6,  1)  =  0.1095781257D  00. 

THEREFORE  THE  6-TH  VARIABLE  LEAVES  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

0,50,  5,11,  7,  1,  9,51,16,15, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XBt  1)  =  X<  0)  =  -0.6225727047D  01 

XBt  2)  =  X(  50)  =  -0 . 9576774908D  00 

XBt  3)  =  X(  5)  =  0.2438492103D  01 

XBt  4)  =  X(  11)  =  0.283B918209D  01 

XBt  5)  =  X<  7)  =  0.4090106653D  00 

XBt  6)  =  X<  1)  =  0.1 09578 1257D  00 

XBt  7)  =  X<  9)  =  0.5393060698D  00 

XBt  8)  =  Xt  51)  =  0.9576774908D  00 

XBt  9)  =  Xt  16)  =  0.1174241027D  00 

XBt  10)  =  Xt  15)  =  0.1597466119D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.95768D  00 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS! 

THE  1-TH  ROW!  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000D  01  0.10000000D  01  -0.10000000P  01 

0.20000000D  01  -0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  2-TH  ROW!  O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-0.10000000D  01  -O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  3-THROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  4-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01  -O.OOOOOOOOD  00 

0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  5-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01 


I 
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-0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  6-TH  ROUS  O.OOOOOOOOD  00  '  O.OOOOOOOOD  00  O.OOOOOOOOD 

-O.IOOOOOOOD  01  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00 

THE  7-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-0.10000000D  01  0.10000000D  01  J.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  8-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.10000000D  01  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  -O.OOOOOOOOD  00 

THE  9-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0 . 1 OOOOOOOD  01 

THE  10-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

PHASE  1  S  ITERATION  2 
THE  REDUCED  COSTS 


O.OOOOOOOOD  00 

00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

•O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

•O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

0.1 OOOOOOOD  01 


THE 

1-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

2-TH 

REDUCED 

COST 

= 

-0, 4430000000D 

01 

THE 

3-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

4-TH 

REDUCED 

COST 

= 

0.4430000000D 

01 

THE 

5-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

6-TH 

REDUCED 

COST 

S 

-O.OOOOOOOOOOD 

00 

THE 

7-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

8-TH 

REDUCED 

COST 

S 

0.1000000000D 

01 

THE 

9-TH 

REDUCED 

COST 

s 

-O.OOOOOOOOOOD 

00 

THE 

10-TH 

REDUCED 

COST 

s 

-O.OOOOOOOOOOD 

00 

THE 

11-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

12-TH 

REDUCED 

COST 

= 

-O.OOOOOOOOOOD 

00 

THE 

J  3-TH 

REDUCED 

COST 

s 

-0.1000000000D 

03 

THE 

14-TH 

REDUCED 

COST 

z 

-O.OOOOOOOOOOD 

00 

THE 

15-TH 

REDUCED 

COST 

r 

-O.OOOOOOOOOOD 

00 

THE 

16-TH 

REDUCED 

COST 

s 

-O.OOOOOOOOOOD 

00 

4)  IS 

THE  NON-BASIC 

VARIABLE  WITH  THE  SMALLEST  INDEX 

x<  _ _ _ _ 

WHICH  HAS  THE  LARGEST  POSITIVE  REDUCED  COST. 
THEREFORE  X<  4)  UILL  ENTER  THE  BASIS. 

THE  COMPONENTS  OF  Y(-»  4) 


Y( 

3, 

4)  = 

-0.4310000000D  01 

Y< 

4, 

4)  = 

0.3480000000D  01 

Y( 

5, 

4)  = 

-0.2410000000D  01 

Y( 

6, 

4)  = 

-0. 5430000000D  01 

Y( 

7, 

4)  = 

0. 1 160000000D  01 

Y( 

8, 

4)  = 

0.4430000000D  01 

Y( 

9, 

4)  = 

0.5430000000D  01 

Y( 

10, 

4)  * 

0.1000000000D  01 

THE  VECTOR  To  LEAVE  THE  BASIS  IS  NOH  DETERMINED. 


7>/Y( 

8) /Y( 

9) /Y( 
10)/Y( 
9)/Y( 


7. 

8, 
9i 

10. 

9. 


4) 

4) 

4) 

4) 

4) 

4) 

4)  = 


0.2162506495D-01. 


THE  XB( I )/Y (If  4)  RATIOS  WITH  Y(I,  4)  >  0.  ARES 
0.81578B  00=  XB<  41/YT  4, 

0.46492P  00=  XB< 

0.21618D  00=  XB( 

0.21625D-01=  XB( 

0.15975D  00=  XB( 

THE  MINIMUM  RATIO  IS  XB< 

THEREFORE  THE  9-TH  VARIABLE  LEAVES  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

0,50,  5,11,  7,  1,  9,51,  4,15, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB(  1)  =  XT  0)  =  -0.62707071S2D  01 

XB<  2)  =  XT  50)  =  -0.86187F.530D  00 

XBT  3)  =  XT  5)  =  0.2531696133D  01 

XBT  4)  =  XT  11)  =  0.27636629B3D  01 

XBT  5)  =  XT  7)  =  0.4611270718D  00 

XBT  6)  =  XT  1)  =  0.2270022283D  00 

XBT  7)  =  XT  9)  =  0.5142209945D  00 

XBT  8)  =  XT  51)  =  0.8M8784530D  00 

XBT  9)  =  XT  4)  =  0.2142506495B-01 

XBT  10)  =  XT  15)  =  0.1381215470D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.86188D  00 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS* 

THE  1-TH  ROW*.  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000D 

0.16169429D  01  -O.IOOOOOOOD  01  -O.OOOOOOOOD  00  0.38305709D  00 

THE  2-TH  ROW!  O.OOOOOOOOD  00  O.IOOOOOOOD  01  -O.OOOOOOOOD 

-0.18416206D  00  -O.OOOOOOOOD  00  O.IOOOOOOOD  01  -0.81583794D  00 

THE  3-TH  ROW.*  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  O.IOOOOOOOD 

-0.20626151D  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.79373849D  00 

THE  4-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.35911602!)  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.64088398D  00 

THE  5-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

— 0.55616943D  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.44383057D  00 

THE  6-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-0.69388939D-14  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01 

THE  7-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-0.12136280D  01  O.IOOOOOOOD  01  -O.OOOOOOOOD  00  0.21362799D  00 

THE  8-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.18416206D  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  0.81583794D  00 


01  O.IOOOOOOOD  01 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
01  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.IOOOOOOOD  01 

•O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 


O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


O.OOOOOOOOD  00 


-23- 

THE  9-TH  ROUT  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

0.18414204D  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0.18414204D  00 

THE  10-TH  ROW*.  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-0.18414204D  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.18414204D  00 

PHASE  1  :  ITERATION  3 
THE  REDUCED  COSTS 


00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

0.10000000D  01 


THE  1-TH  REDUCED  COST 
THE  2-TH  REDUCED  COST 
THE  3-TH  REDUCED  COST 
THE  4-TH  REDUCED  COST 
THE  5-TH  REDUCED  COST 
THE  4-TH  REDUCED  COST 
THE  7-TH  REDUCED  COST 
THE  8-TH  REDUCED  COST 
THE  9-TH  REDUCED  COST 
THE  10-TH  REDUCED  COST 
THE  11-TH  REDUCED  COST 
THE  12-TH  REDUCED  COST 
THE  13-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
THE  15-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
X<  8)  IS  THE  NON-BASIC  VARIABLE  WITH  THE  SMALLEST 
WHICH  HAS  THE  LARGEST  POSITIVE  REDUCED  COST. 
THEREFORE  X(  8)  WILL  ENTER  THE  BASIS. 

THE  COMPONENTS  OF  Y<-,  8) 


-O.OOOOOOOOOOD  00 
-0.2220444049D-15 
-O.OOOOOOOOOOD  00 
0.2220444049D-15 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
0.1841420424D  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-0.1841420424D  00 
-O.OOOOOOOOOOD  00 
-0.8158379374D  00 
-O.OOOOOOOOOOD  00 


INDEX 


Y( 

3, 

8)  = 

-0.2042415101D  OO 

Y( 

4, 

8)  = 

0.3591 140221D  00 

Y( 

5, 

8)  = 

-0.5541694291B  00 

Y( 

4, 

8)  = 

-0.4938893904D-14 

Y( 

7, 

8)  = 

-0.1213427993B  01 

Y( 

8, 

8)  = 

0.1841420424D  00 

Y< 

9, 

8)  = 

0.1841420424D  00 

Y( 

10, 

8)  = 

-0.1841420424D  00 

THE  VECTOR  TO  LEAVE  THE  BASIS  IS  NOW  DETERMINED. 


8)  > 
8) 
8) 
8) 

8)  * 


o.  are: 


0.1174241027B  00. 


THE  XB(I)/Y(Ii  8)  RATIOS  WITH  Y(I, 

0.74957D  01=  XB<  4>/Y<  4, 

0.44800D  01=  XB<  8)/Y(  8, 

0.11742D  00=  XB(  9>/Y<  9, 

THE  HINIhUM  RATIO  IS  XB<  9)/Y<  9, 

THEREFORE  THE  9-TH  VARIABLE  LEAVES  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

0,50,  5,11,  7,  1,  9,51,  8,15, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB(  1)  =  X<  0)  =  -0.4577999355D  01 

XB(  2)  =  X(  50)  =  -0.8402533881D  00 

XB(  3)  =  X(  5)  =  0. 25559 14205D  01 

XB(  4)  =  X<  11)  =  0.2721494107D  01 

XB<  5)  =  X(  7)  =  0.5244347480D  00 

XB(  4)  =  X(  1)  =  0.22700222B3D  00 

XB(  7)  =  X<  9)  =  0. 4547301 725D  00 

XB<  8)  =  X<  51)  =  0.8402533881D  00 

XB<  9)  =  X(  8)  =  0.1174241027D  00 

XB(  10)  =  X(  15)  =  0.1597444U9D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.84025D  00 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS: 

THE  1-TH  ROW:  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000D 

-0.10000000D  01  -0.10000000D  01  -O.OOOOOOOOD  00  0.30000000D  01 

THE  2-TH  ROW!  O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01  -0.10000000D  01 

THE  3-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01 

THE  4-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01 

THE  5-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01 

THE  6-TH  R0W5  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01 

THE  7-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000D  01 

THE  8-TH  ROWS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  0.10000000D  01 

THE  9-TH  ROWS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

O.IOOOOOOOD  01  O.OOOOOOOOD  00  O.OOOOOOOOD  00  -O.IOOOOOOOD  01 

THE  10-TH  R8W'  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

PHASE  1  5  ITERATION  4 
THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  2-TH  REDUCED  COST  =  0.1000000000D  01 

THE  3-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 


01  O.IOOOOOOOD  01 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
01  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.IOOOOOOOD  01 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.IOOOOOOOD  01 


-O.OOOOOOOOD  00 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


4- TH  REDUCED  COST 

5- TH  REDUCED  COST 
4-TH  REDUCED  COST 

7- TH  REDUCED  COST 

8- TH  REDUCED  COST 

9- TH  RETHJCED  COST 

10- TH  REDUCED  COST 

11- TH  REDUCED  COST 

12- TH  REDUCED  COST 

13- TH  REDUCED  COST 

14- TH  REDUCED  COST 

15- TH  REDUCED  COST 

16- TH  REDUCED  COST 
X(  2)  IS  THE  NON-BASIC  VARIABLE  WITH  THE  SMALLEST  INDEX 
WHICH  HAS  THE  LARGEST  POSITIVE  REDUCED  COST. 

THEREFORE  X<  2)  WILL  ENTER  THE  BASIS. 

THE  COMPONENTS  Of  Y(-,  2) 


THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


-O.IOOOOOOOOOD  01 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.IOOOOOOOOOD  01 
-O.OOOOOOOOOOD  00 


Y( 

3, 

2)  = 

-0.1 120000000D  01 

Y< 

4, 

2)  = 

0.1950000000D  01 

Y( 

5, 

2)  = 

-0.3020000000D  01 

Y( 

4, 

2)  = 

O.OOOOOOOOOOD  00 

Y( 

7, 

2)  = 

-0.6590000000D  01 

Y< 

8, 

2)  = 

O.IOOOOOOOOOD  01 

Y( 

9, 

2)  = 

-0.5430000000D  01 

Y( 

10, 

2)  = 

-O.IOOOOOOOOOD  01 

THE  VECTOR  TO  LEAVE  THE  BASIS  IS  NOW  DETERMINED. 

THE  XB(I)/Y<I,  2)  RATIOS  WITH  YII,  2)  >  0.  ARE? 

0.13956D  01=  XB(  4>/Y(  4.  2) 

0.84025D  00=  XB(  8>/Y(  8,  2) 

THE  MINIMUM  RATIO  IS  XB<  8)/Y(  8.  2)  =  0.8402533881B  00. 

THEREFORE  THE  8-TH  VARIABLE  !  PAVES  THE  BASIS. 

THE  BASIC  VARIABLES  ARE  NOW 

0.50,  5,11,  7,  1,  9,  2,  8,15, 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 


XB( 

XB( 

XB( 

XB< 

XB< 

XB( 

XB( 

XB< 

XB( 

XB( 


1) 

2) 

3) 

4) 

5) 

6) 

7) 

8) 

9) 

10) 


X(  0) 
X(  50) 


X< 

X( 

X( 

X( 

X< 

X( 

X( 

X( 


5) 

11) 

7) 
1) 
9) 
2) 

8) 
15) 


-O.OOOOOD  00 


-0.1851800000D  02 
-O.OOOOOOOOOOD  00 
0.3497000000D  01 
0.1083000000D  01 
0.3064000000D  01 
0.2270022283D  00 
0.4194000000D  01 
0.84 0253388 ID  00 
0.4680000000D  01 
O.IOOOOOOOOOD  01 
THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS: 

THE  1-TH  ROW!  0.10000000D  01  -O.OOOOOOOOB  00  -0.10000000D 

0.10000000D  01  -0.10000000D  01  0.14210000D  02  -0.11210000D  02 

THE  2-THROW.*  O.OOOOOOOOB  00  0.10000000D  01  -O.OOOOOOOOB 

O.OOOOOOOOB  00  -O.OOOOOOOOB  00  -O.OOOOOOOOB  00  -O.OOOOOOOOB  00 

THE  3-THROW:  O.OOOOOOOOB  00  -O.OOOOOOOOB  00  0.10000000D 

O.OOOOOOOOB  00  -O.OOOOOOOOB  00  -0.11200G00D  01  0.12000000D  00 

THE  4-TH  ROW!  O.OOOOOOOOB  00  -O.OOOOOOOOB  00  -O.OOOOOOOOB 

O.OOOOOOOOB  00  -O.OOOOOOOOB  00  0.19500000D  01  -0.95000000D  00 

THE  5-TH  ROW!  O.OOOOOOOOB  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

O.OOOOOOOOB  00  -O.OOOOOOOOD  00  -0.30200000D  01  0.20200000D  01 

THE  4-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01 

THE  7-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

O.OOOOOOOOD  00  0.10000000D  01  -0.45900000D  01  0.55900000D  01 

THE  8-TH  ROW!  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0.10000000D  01  0.100000000  01 

THE  9-TH  ROWS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.10000000D  01  -O.OOOOOOOOD  00  -0.54300000D  01  0.44300000D  01 

THE  10-TH  ROW.'  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01  0.10000000D  01 

PHASE  1  S  ITERATION  5 
THE  REDUCED  COSTS 


01  0.10000000D  01 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
01  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
00  -0.10000000D  01 

-O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

0.10000000D  01 


THE 

THE 

THE 

TIC 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


1- TH  REDUCED  COST  = 

2- TH  REDUCED  COST  = 

3- TH  REDUCED  COST  = 

4- TH  REDUCED  COST  = 

5- TH  REDUCED  COST  = 
4-TH  REDUCED  COST  = 

7- TH  REDUCED  COST  = 

8- TH  REDUCED  COST  = 

9- TH  REDUCED  COST  = 

10- TH  REDUCED  COST  = 

11- TH  REDUCED  COST  = 

12- TH  REDUCED  COST  = 

13- TH  REDUCED  COST  = 


-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 
-O.OOOOOOOOOOD  00 


-0. 1 OOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
0.1 OOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


THE  14-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  15-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  VALUE  OF  THE  SUM  OF  THE  ARTIFICIAL  VARIABLES  AT  THE  EK'D  OF  PHASE  1 

IS  ZERO.  (ACTUALLY  IT  IS  O.OOOOOOOOOOD  00) 

PHASE  2  NOV  BEGINS. 


THE  INITIAL  FEASIBLE  SOLUTION  IS  AS  FOLLOWS! 

X<  1)  =  0.2270022283D  00 

X<  2)  =  0.84025338811*  00 

X<  3)  =  O.OOOOOOOOOOD  00 

X<  4)  =  O.OOOOOOOOOOD  00 

X(  5)  =  0 . 3497000000D  01 

X(  6)  =  O.OOOOOOOOOOD  00 

X<  7)  =  0.3064000000D  01 

X(  8)  =  0.4680000000D  01 

X(  9)  =  0.6194000000D  01 

X(  10)  =  O.OOOOOOOOOOD  00 

X(  11)  =  0.1083000000D  01 

X(  12)  =  O.OOOOOOOOOOD  00 

X(  13)  =  O.OOOOOOOOOOD  00 

X(  14)  =  O.OOOOOOOOOOD  00 

X(  15)  =  0.1000000000D  01 

X(  16)  =  O.OOOOOOOOOOD  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.18518D  02 

PHASE  2  1  ITERATION  1 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.2220446049B-15 

THE  2-TH  REDIKCED  COST  =  -O.OOOOOOOOOOD  00 

THE  3-TH  REDUCED  COST  =  0.2220446049D-15 

THE  4-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  5-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  6-TH  REDUCED  COST  =  -0.2000000000D  01 

THE  7-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  10-TH  REDUCED  COST  =  -0.2000000000D  01 

THE  1 1-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

THE  12-TH  REDUCED  COST  =  -0.2000000000D  01 

THE  13-TH  REDUCED  COST  =  -0.2000000000D  01 

THE  14-TH  REDUCED  COST  =  -0.2000000000D  01 

THE  15-TH  REDUCED  COST  =  -0.11210000001*  02 

THE  16-TH  REDUCED  COST  =  -O.OOOOOOOOOOD  00 

ALL  OF  THE  REDUCED  COSTS  ARE  NON-POSITIVE. 

THEREFORE  THE  CURRENT  BASIC  FEASIBLE  SOLUTION  IS  AN  OPTIMAL  SOLUTION. 

THE  OPTIMAL  SOLUTION  IS  AS  FOLLOWS! 

X<  1)  =  0.2270022283D  00 

X<  2)  =  0.8402533881D  00 

X<  3)  =  O.OOOOOOOOOOD  00 

X<  4)  *  O.OOOOOOOOOOD  00 

X<  5)  =  0.3497000000D  01 

X(  6)  =  O.OOOOOOOOOOD  00 

X<  7)  =  0.3064000000D  01 

X(  8)  *  0.4680000000D  01 

X(  9)  =  0.6194000000D  01 

X(  10)  =  O.OOOOOOOOOOD  00 

X(  11)  =  0. 10830000001*  01 

X(  12)  *  O.OOOOOOOOOOD  00 

X<  13)  =  O.OOOOOOOOOOD  00 

X(  14)  =  O.OOOOOOOOOOD  00 

X(  15)  =  0.1000000000D  01 

X(  16)  *  O.OOOOOOOOOOD  00 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.18518B  02 


THE  AUXILIARY  LEAST  SQUARES  ESTIMATE.  BETAO.  OF  THE  REGRESSION  PARAMETER  VECTOR.  BETA 
LEAST  SQUARES  ESTIMATE  OF  BETA(  1)  =  3.227002 

LEAST  SQUARES  ESTIMATE  OF  BETA!  2)  =  -1.159747 

MR.A'S  ANSWER  1  THE  ESTIMATE  OF  THE  REGRESSION  PARAMETER  VECTOR  WHICH  MINIMIZES  THE  SUH  OF  THE  ABSOLUTE  RESIDUALS 

LI  ESTIMATE  OF  BETA<  1)  =  3.000000 

LI  ESTIMATE  OF  BETA!  2)  «  -2.000000 

THE  RESIDUALS  1 

THE  1-TH  RESIDUAL  =  3.497000 

THE  2-TH  RESIDUAL  *  -1.083000 

THE  3-TH  RESIDUAL  =  3.064000 

THE  4-TH  RESIDUAL  >  4.680000 

THE  5-TH  RESIDUAL  =  6.194000 

THE  SUM  OF  THE  ABSOLUTE  RESIDUALS  =  18.518000 

THE  MAXIMUM  ABSOLUTE  RESIDUAL  *  6.194000 


THE  MONTE  CARLO  STUDY  TO  ESTIMATE  THE  COVARIANCE  OF  THE  LI  ESTIMATE  OF  BETA  NOW  BEGINS 
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RESULTS  FOR  MONTE  CARLO  SAMPLE  NUMBER 
THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0,810298  -0.168497 

-0.168497  0.046521 

THE  VALUE  OF  BET AO  TO  COMPUTE  RHS 
-0.594250  0.136933 

THE  RIGHT  HAND  SIDES  FOR  THIS  SAMPLE 


i: 


BRHS1(  1) 


0.000000 


BRHSK  2)  = 

0.000000 

BRHSK  3)  = 

0.067696 

BRHSK  4)  = 

-0.244174 

BRHSK  5)  = 

0.226595 

BRHSK  6)  = 

-0.039891 

BRHSK  7)  = 

-0.010226 

BRHSK  8)  = 

0.236092 

BRHSK  9)  = 

0.582285 

BRHSK  10)  = 

-0.184188 

THE  LEAST  SQUARES  ESTIMATE  OF  BETA  IN  THIS 

BETA!  1)  = 

-0.594250 

BETA<  2)  = 

0.136933 

THE  EPSILONS  GENERATED  FOR  THIS  SAMPLE 

EPSILON!  1)  = 

-0.373189 

EPSILON!  2)  = 

-0.571406 

EPSILON!  3)  = 

0.045882 

EPSILON!  4)  = 

0.109403 

EPSILON!  5)  = 

0.297911 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 

XB<  0)  = 

-0.36609D  01 

XB!  1)  = 

O.OOOOOD  00 

XB<  2)  = 

-0.12685D  00 

XB!  3)  = 

0.15138D  00 

XB!  4)  = 

0.68981D  00 

XB!  5)  = 

-0.58228D  00 

XB!  6)  = 

0.16889D  01 

XB!  7)  = 

0.34619D  00 

XB!  8)  = 

0.12577D  01 

XB!  9)  = 

0.16200D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 
THE  REDUCED  COSTS 


-0.36609D  01 


THE  1-TH  REDUCED  COST 
THE  2-TH  REDUCED  COST 
THE  3-TH  REDUCED  COST 
THE  4-TH  REDUCED  COST 
THE  5-TH  REDUCED  COST 
THE  6-TH  REDUCED  COST 
THE  7-TH  REDUCED  COST 
THE  8-TH  REDUCED  COST 
THE  9-TH  REDUCED  COST 
THE  10-TH  REDUCED  COST 
THE  11-TH  REDUCED  COST 
THE  12-TH  REDUCED  COST 
THE  13-TH  REDUCED  COST 
THE  14-TH  REDUCED  COST 
THE  15-TH  REDUCED  COST 
16-TH  REDUCED  COST 


-O.OOOOOD  00 
-O.OOOOOD  00 
0.22204D-15 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.20000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.20000D  01 
-O.OOOOOD  00 
-0.20000D  01 
-0.20000D  01 
-0.20000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 


THE 

TENTATIVELY  THE 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


6-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


IT  CAN  PE  ZERO  OTHERWISE  TOO,) 

YR( 

1) 

S 

O.OOOOOD  00 

YR( 

2) 

= 

O.OOOOOD  00 

YR< 

3) 

s 

-0.10000D  01 

YR( 

4) 

s 

O.OOOOOD  00 

YR< 

5) 

s 

O.OOOOOD  00 

YR! 

6) 

s 

O.OOOOOD  00 

YR< 

7) 

« 

O.OOOOOD  00 

YR( 

8) 

s 

O.OOOOOD  00 

YR( 

9) 

s 

O.OOOOOD  00 

YR! 

10) 

s 

O.OOOOOD  00 

YR! 

11) 

r 

O.OOOOOD  00 

YR! 

12) 

s 

O.OOOOOD  00 

YR! 

13) 

z 

O.OOOOOD  00 

YR! 

14) 

s 

O.OOOOOD  00 

YR! 

15) 

s 

O.OOOOOD  00 

THE  RATIO1^ 

THE 

■MH  ft 

THE 


_  _ .........  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  Tiff  3-TH  NET  PRICE  TO  YR<  3)  = 
6-TH  VARIABLE  IS  LEAVING  THE  BASIS, 


0.22204D-15 
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THE  3-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOUSS 
THE  1-TH  ROUT  O.IOOOOOOOD  01  -O.OOOOOOOOD  00  -0.10000000D 

-0.10000000D  01  -O.IOOOOOOOD  01  0.H210000D  02  -0.11210000D  02 

THE  2-TH  ROW!  O.OOOOOOOOD  00  O.IOOOOOOOD  01  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  3-TH  ROU!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  O.IOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.11200000B  01  0.12000000D  00 

THE  4-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.19500000D  01  -0.95000000D  00 

THE  5-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.30200000D  01  0.20200000D  01 

THE  4-TH  ROUS  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.IOOOOOOOD  01 

THE  7-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  O.IOOOOOOOD  01  -0.65900000D  01  0.55900000D  01 

THE  8-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01 

THE  9-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

O.IOOOOOOOD  01  -O.OOOOOOOOD  00  -0,543000000  01  0.44300000D  01 

THE  10-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01 

THE  BASIC  VARIABLES  ARE  NOU 

0,  50.  5.  11»  7»  3*  9*  2.  Si  15. 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOU 


01  O.IOOOOOOOD  01 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

•O.OOOOOOOOD  00 
01  -O.OOOOOOOOD  00 

•O.OOOOOOOOD  00 
00  -O.IOOOOOOOD  01 

•O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.IOOOOOOOD  01 


XBI 

1)  = 

-O.OOOOOD  00 

XBI 

2)  = 

-0.12685D  00 

XBI 

3)  = 

0. 151380  00 

XBI 

4)  = 

0.68981D  00 

XBI 

5)  = 

0.58228D  00 

XBI 

6)  = 

0.16889D  01 

XBI 

7)  = 

0.34419D  00 

XBI 

8)  = 

0.12577D  01 

XBI 

9)  = 

0.16200D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 
THE  REDUCED  COSTS 


-0.34409D  01 


THE 

1-TH 

REDUCED 

COST 

r 

-O.OOOOOD 

00 

THE 

2-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

THE 

3-TH 

REDUCED 

COST 

r 

-O.OOOOOD 

00 

THE 

4-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

THE 

5-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

THE 

4-TH 

REDUCED 

COST 

r 

-0.20000D 

01 

THE 

7-TH 

REDUCED 

COST 

r 

-O.OOOOOD 

00 

THE 

8-TH 

REDUCED 

COST 

r 

-O.OOOOOD 

00 

THE 

9-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

THE 

10-TH 

REDUCED 

COST 

r 

-0.20000D 

01 

THE 

11— TH 

REDUCED 

COST 

s 

-O.OOOOOD 

00 

THE 

12-TH 

REDUCED 

COST 

= 

-0.20000D 

01 

THE 

13-TH 

REDUCED 

COST 

= 

-0.20000P 

01 

THE 

14-TH 

REDUCED 

COST 

= 

-0.20000D 

01 

THE 

15-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

THE 

14-TH 

REDUCED 

COST 

= 

-O.OOOOOD 

00 

TENTATIVELY  THE 

3-TH  BASIC 

VARIABLE  IS  LEAVING 

THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 

IT  CAN  BE  ZERO  OTHERUISE  TOO.) 

YR(  1)  =  O.OOOOOD  00 

YR(  2)  =  O.OOOOOD  00 

YR(  3)  =  O.OOOOOD  00 

YR(  4)  =  O.OOOOOD  00 

YR(  5)  =  O.OOOOOD  00 

YR(  6)  =  O.OOOOOD  00 

YR(  7)  =  O.OOOOOD  00 

YR<  B)  =  O.OOOOOD  00 

YR(  9)  =  O.OOOOOD  00 

YRI  10)  =  -0.10000D  01 

YR(  11)  =  O.OOOOOD  00 

YR(  12)  =  O.OOOOOD  00 

YR<  13)  *  O.OOOOOD  00 

O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  10-TH  NET  PRICE  TO  YR(  10)  =  -0.20000D  01 

THE  3-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  10-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOUSS 
THE  1-TH  ROUS  O.IOOOOOOOD  01  -O.OOOOOOOOD  00  O.IOOOOOOOD 

-O.IOOOOOOOD  01  -O.IOOOOOOOD  01  0.11970000D  02  -0.10970000D  02 

THE  2-TH  ROUS  O.OOOOOOOOD  00  O.IOOOOOOOD  01  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 


YR(  14) 
YR<  15) 
YR{  14)  = 


01  O.IOOOOOOOD  01 
O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.IOOOOOOOD  01 
O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 


THE  3-TH  ROW!  O.DOOOOOOOD  00  O.OOOOOOOOD  00  -0. 100000000 

0 >000000000  00  O.OOOOOOOOD  00  0.11200000D  01  -0.120000000  00 

THE  4-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.19500000D  01  -0.95000000D  00 

THE  5-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0. 30200000D  01  0.20200000D  01 

THE  3-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01 

THE  7-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  0.10000000D  01  -0.&5900000D  01  0.55900000D  01 

THE  8-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01  0.10000000D  01 

THE  9-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

0.10000000D  01  -O.OOOOOOOOD  00  -0.54300000D  01  0.44300000D  01 

THE  10-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.10000000D  01  0.10000000D  01 

THE  6ASIC  VARIABLES  ARE  NOW 

Oi  50t  lOi  lit  7t  3t  9t  2t  8t  15t 
THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOU 
XB(  1)  =  -O.OOOOOD  00 

XB(  2)  =  0.12685D  00 

XB(  3)  =  0.15138D  00 

XB(  4)  =  0.6B981D  00 

XB(  5)  =  0.58228D  00 

XB(  6)  =  0.16B89D  01 

XB(  7)  =  0.34&19D  00 

XB(  8)  =  0.12577D  01 

XB(  9)  =  0.16200D  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.39146D  01 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -0.20000D  01 

THE  3-TH  REDUCED  COST  =  -0.20000D  01 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  12-TH  REDUCED  COST  =  -0.20000D  01 

THE  13-TH  REDUCED  COST  =  -0.20000D  01 

THE  14-TH  REDUCED  COST  =  -0.20000D  01 

THE  15-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTINAL. 

THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS: 

X<  0)  =  -0.39146D  01 

X<  50)  =  -O.OOOOOD  00 

X<  10)  *  0.12685D  00 

X(  11)  =  0.15138D  00 

X(  7)  =  0.68981D  00 

X(  3)  =  0.58228D  00 

X(  9)  =  0.16889D  01 

X(  2)  =  0.34619D  00 

X<  8)  =  0.12577D  01 

X(  15)  =  0.13200D  00 

THE  OPTINAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.39143D  01 

RESULTS  FOR  MONTE  CARLO  SAMPLE  NUMBER  21 
THE  PROBLEM  PI  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.810298  -0.138497 

-0.138497  0.043521 

THE  VALUE  OF  BET AO  TO  COMPUTE  RHS 
0.217174  0.114131 

THE  RIGHT  HAND  SIDES  FOR  THIS  SAMPLE 
8RHSH  1)  =  0.000000 

BRHSK  2)  =  0.000000 

BRHSH  3)  =  0.534109 

BRHSK  4)  *  0.157815 

BRHSK  5)  *  -0.919343 

BRHSK  3)  *  -0.320504 

BRHSK  7)  *  0.547923 

BRHSK  8)  «  0.070373 

BRHSK  9)  *  0.294137 

BRHSK  10)  »  -0.117734 

THE  LEAST  SQUARES  ESTIMATE  OF  BETA  IN  THIS  SAMPLE 
BETA!  1)  *  0.217174 

BETA!  2)  *  0.114131 


01  O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -0.10000000D  01 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

•O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

■O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

O.OOOOOOOOD  00 
00  -O.OOOOOOOOD  00 

0.10000000D  01 


O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
0.10000000D  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 


(E 
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BRHSK  3)  =  -1.062358 

BRUSH  0)  =  1.453735 

BRHSK  5)  *  -0.206058 

BRHSK  6)  =  -0.171215 

BRHSK  7)  =  -0.014105 

BRHSK  8)  =  0.358040 

BRHSK  9)  =  0.985521 

BRHSK  10)  *  -0.297676 

THE  LEAST  SQUARES  ESTIMATE  Of  BETA  IN  THIS  SAMPLE 
BETA!  1)  =  0.011806 

BETA!  2)  =  -0.097417 

THE  EPSILONS  GENERATED  FOR  THIS  SAHPLE 
EPSILON!  1)  =  -1.159658 

EPSILON!  2)  =  1.275578 

EPSILON!  3)  =  -0.488451 

EPSILON!  4)  =  -0.688384 

EPSILON!  5)  =  -0.644277 

THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 


(S 


XB( 

0)  = 

-0.91263D  01 

XB< 

1)  = 

O.OOOOOD  00 

XB! 

2)  = 

0.13451D  01 

XB( 

3)  = 

0.16918D  01 

XB! 

4)  = 

0.70341D  00 

XB! 

5)  = 

0.98552D  00 

XB! 

6)  = 

0.31355D  01 

XB! 

7)  = 

0.62748D  00 

XB! 

8)  = 

0.22505D  01 

XB! 

9)  = 

0.32980D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  = 
THE  REDUCED  COSTS 


-0.91263D  01 


1- TH  REDUCED  COST  = 

2- TH  REDUCED  COST  = 

3- TH  REDUCED  COST  = 

4- TH  REDUCED  COST  = 

5- TH  REDUCED  COST  = 

6- TH  REDUCED  COST  = 

7- TH  REDUCED  COST  = 

8- TH  REDUCED  COST  = 

9- TH  REDUCED  COST  = 

10-  TN  REDUCED  COST  = 

11- TH  REDUCED  COST  * 

12- TH  REDUCED  COST  = 

13-  TH  REDUCED  COST  = 

14- TH  REDUCED  COST  = 

15- TH  REDUCED  COST  = 

16- TH  REDUCED  COST  = 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OPTIMAL, 
THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS! 


THE 

THE 

THE 

THE 

T* 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 

THE 


-O.OOOOOD  00 
0.28866D-14 
-0.28866D-14 
-O.OOOOOD  00 
-0.20000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.20000D  01 
-0.20000D  01 
-0.20000D  01 
-0.20000B  01 
-O.OOOOOD  00 
-O.OOOOOD  00 


X! 

0)  = 

-0.91263D  01 

X!  50)  = 

O.OOOOOD  00 

X! 

10)  = 

0.13451D  01 

X! 

6)  = 

0.16918D  01 

X! 

7)  = 

0.70341D  00 

X! 

1)  = 

0.98552D  00 

X! 

9)  = 

0.31355D  01 

X! 

4)  = 

0.62748D  00 

X! 

8)  = 

0.22505D  01 

X!  15)  = 

0.32980D  00 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.91263D  01 

RESULTS  FOR  MONTE  CARLO  SAMPLE  NUMBER  20! 

THE  PROBLEM  P2  HAS  BEEN  SELECTED 
THE  INVERSE  OF  XTX 

0.810298  -0.168497 

-0.168497  0.046521 

THE  VALUE  OF  BETAO  TO  COMPUTE  RHS 
0.173768  -0,207514 

THE  RIGHT  HAND  SIDES  FOR  THIS  SAMPLE 
BRHSK  1)  «  0.000000 

BRHSK  2)  -  0.000000 

BRHSK  3)  *  0.478624 

BRHSK  4)  «  -0.259669 

BRHSK  5)  =  -0.702085 

BRHSK  6)  =  0.942444 

BRHSK  7)  •  -0.459314 

BRHSK  8)  «  0.349361 

BRHSK  9)  *  0.955143 

BRHSK  10)  -  -0.267826 

THE  LEAST  SQUARES  ESTIMATE  OF  BETA  IN  THIS  SAHPLE 
BETA!  1)  *  0.173768 

BETA!  2)  *  -0.207514 


ECO*J  T/AJU  E  O 
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THE  EPSILONS  GENERATED  FOR  THIS  SAMPLE 


EPSILONt 

EPSILONt 

EPSILOHt 

EPSILOHt 

EPSILOHt 


0.419977 

-0.490553 

-1.155009 

-0.010589 

-1.653063 


THE  INITIAL  VALUES  OF  THE  BASIC  VARIABLES 
XBt  0)  =  -0.57911D  01 

XBt  1)  =  O.OOOOOD  00 

XBt  2)  =  -0.20196D  00 

XBt  3)  =  -0.33537D-01 

XBt  4)  =  0.17223D  00 

XBt  5)  =  0.95514D  00 

XBt  6)  =  0.25776D  01 

XBt  7)  =  0.60578D  00 

XBt  8)  =  0.32767D  01 

XBt  9)  =  0.33796D  00 

THE  INITIAL  VALUE  OF  THE  OBJECTIVE  FUNCTION 
THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.28866D 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD 

THE  4-TH  REDUCED  COST  =  0.28866D 

THE  5-TH  REDUCED  COST  =  -0.20000D 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD 

TK  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  10-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  11-TH  REDUCED  COST  =  -0.20000D  01 

THE  12-TH  REDUCED  COST  =  -0.20000D  01 

THE  13— TH  REDUCED  COST  =  -0.200000  01 

THE  14-TH  REDUCED  COST  =  -0.20000D  01 

THE  15-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

TENTATIVELY  THE  3-TH  BASIC  VARIABLE  IS  LEAVING  TH 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  Tiff  BASIS. 


-0.57911D  01 


-0.28866D-14 
-O.OOOOOD  00 
-O.OOOOOD  00 
0.28866D-14 
-0.20000D  01 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-O.OOOOOD  00 
-0.20000D  01 
-0.20000D  01 
-0.200000  01 
-0.200000  01 
-O.OOOOOD  00 
-O.OOOOOD  00 


3-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

YRt 

1) 

r 

O.OOOOOD  00 

YR( 

2) 

s 

O.OOOOOD  00 

YRt 

3) 

s 

O.OOOOOD  00 

YRt 

4) 

= 

O.OOOOOD  00 

YRt 

5) 

s 

-0.10000D  01 

YRt 

6) 

2 

O.OOOOOD  00 

YRt 

7) 

2 

O.OOOOOD  00 

YRt 

8) 

2 

O.OOOOOD  00 

YRt 

9) 

2 

O.OOOOOD  00 

YRt 

10) 

2 

O.OOOOOD  00 

YRt 

11) 

2 

O.OOOOOD  00 

YRt  12) 

2 

O.OOOOOD  00 

YRt 

13) 

2 

O.OOOOOD  00 

YRt 

14) 

2 

O.OOOOOD  00 

YRt  15) 

2 

O.OOOOOD  00 

YRt 

16) 

2 

O.OOOOOD  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  5-TH  NET  PRICE  TO  YR<  5)  =  -0.20000D  01 

THE  3-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  5-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS? 

THE  1-TH  ROW!  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000D  01  -0.10000000D  01 

-0.10000000D  01  -0.10000000D  01  0.18110000D  02  -0.13110000D  02  -O.OOOOOOOOD  00 

THE  2-TH  RON*  O.OOOOOOOOD  00  O.IOOOOOOOD  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  3-TH  ROW:  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.IOOOOOOOD  01  O.OOOOOOOOD  00 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0.1 1200000D  01  0.12000000D  00  O.OOOOOOOOD  00 

THE  4-TH  R0U1  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  O.IOOOOOOOD  01 

-0.000000000  00  -O.OOOOOOOOD  00  -0.19500000D  01  0.95000000D  00  -O.OOOOOOOOD  00 

THE  5-TH  RtMS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.30200000D  01  0.20200000D  01  -O.OOOOOOOOD  00 

THE  6-TH  ROUS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  O.IOOOOOOOD  01  -O.OOOOOOOOD  00 

THE  7-TH  ROD:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  'O.IOOOOOOOD  01  -0.6S900000D  01  0.55900000D  01  -O.OOOOOOOOD  00 

THE  8-TH  ROVJ  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01  -O.OOOOOOOOD  00 

THE  9-TH  RDM:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

O.IOOOOOOOD  01  -O.OOOOOOOOD  00  -0.54300000D  01  0.44300000D  01  -O.OOOOOOOOD  00 

THE  10-TH  ROW!  0.00000000!  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01  O.IOOOOOOOD  01 

Tiff  BASIC  VARIABLES  ARE  NOV 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
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0*  50*  5«  6*  7*  3»  9»  2*  8>  IS* 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB<  1)  =  -O.OOOOOD  00 

XB<  2)  *  0.2019&D  00 

XB(  3)  =  -0.33537D-01 

XB<  4)  *  0.17223D  00 

XB(  5)  *  0.95514D  00 

XB(  6)  =  0.2577AD  01 

XB<  7)  *  0.60578D  00 

XB(  8)  =  0.32767D  01 

XB(  9>  =  0 >337960  00 

THE  CURRENT  VALUE  OF  THE  OBJECTIVE  FUNCTION  =  -0.619500  01 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  =  -0.28866D-14 

THE  2-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  6-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  8-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  9-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  10-TH  REDUCED  COST  =  -0.20000D  01 

THE  11-TH  REDUCED  COST  =  -0.20000D  01 

THE  12-TH  REDUCED  COST  =  -0.20000D  01 

THE  13-TH  REDUCED  COST  =  -0.20000D  01 

THE  H-TH  REDUCED  COST  =  -0.20000D  01 

THE  15-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  16-TH  REDUCED  COST  =  -O.OOOOOD  00 

TENTATIVELY  THE  4-TH  BASIC  VARIABLE  IS  LEAVING  THE  BASIS 
THE  YRJ'S 

(  YR(J)  =  0  IF  THE  J-TH  VARIABLE  IS  IN  THE  BASIS. 


IT  CAN  BE  ZERO  OTHERWISE  TOO.) 

O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
-0.10000D  01 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 
O.OOOOOD  00 

THE  RATIO  OF  THE  NET  PRICES  TO  THE  NEGATIVE  YRJ'S 

THE  RATIO  OF  THE  11-TH  NET  PRICE  TO  YR(  11)  =  -0.20000D  01 

THE  4-TH  VARIABLE  IS  LEAVING  THE  BASIS. 

THE  11-TH  VARIABLE  IS  ENTERING  THE  BASIS. 

THE  CURRENT  BASIS  INVERSE  IS  AS  FOLLOWS. 

THE  1-TH  ROW!  0.10000000D  01  -O.OOOOOOOOD  00  -0.10000000B  01  0.10000000D  01 

-0.10000000D  01  -0.10000000D  01  0.14210000D  02  -0.11210000D  02  -O.OOOOOOOOD  00 

THE  2-TH  ROW!  O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

THE  3-TH  ROW?  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.11200000D  01  0.12000000D  00  -O.OOOOOOOOD  00 

THE  4-TH  ROW!  O.OOOOOOOOD  00  O.OOOOOOOOD  00  O.OOOOOOOOD  00  -0.10000000D  01 

O.OOOOOOOOD  00  O.OOOOOOOOD  00  0.19500000D  01  -0.95000000D  00  O.OOOOOOOOD  00 

THE  5-TH  ROW;  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -0.30200000D  01  0.20200000D  01  -O.OOOOOOOOD  00 

THE  6-TH  ROWS  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  0.10000000D  01  -O.OOOOOOOOD  00 

THE  7-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  0.10000000D  01  -0.65900000D  01  0.55900000D  01  -O.OOOOOOOOD  00 

THE  8-TH  ROW!  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01  -O.OOOOOOOOD  00 

THE  9-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

O.IOOOOOOOD  01  -O.OOOOOOOOD  00  -0.54300000D  01  0.44300000D  01  -O.OOOOOOOOD  00 

THE  10-TH  ROW:  O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.OOOOOOOOD  00 

-O.OOOOOOOOD  00  -O.OOOOOOOOD  00  -O.IOOOOOOOD  01  O.IOOOOOOOD  01  O.IOOOOOOOD  01 

THE  BASIC  VARIABLES  ARE  NOW 

0'  50*  5*  11*  7*  3*  9*  2*  8*  15* 

THE  VALUES  OF  THE  BASIC  VARIABLES  ARE  NOW 
XB<  1)  *  -O.OOOOOD  00 


YR( 

1)  = 

YR< 

2)  * 

YR( 

3)  = 

YR< 

4)  = 

YR( 

5)  = 

YR< 

6)  = 

YR( 

7)  = 

YR( 

8)  = 

YR( 

9)  * 

YR( 

10)  = 

YR( 

11)  = 

YR( 

12)  = 

YR< 

13)  = 

YR< 

14)  = 

YR( 

15)  = 

YR< 

16)  = 

XB<  2)  *  0.20196D  00 


XB(  3)  *  0.335370-01 


XB<  4)  *  0.17223D  00 


-O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
O.OOOOOOOOD  00 
O.IOOOOOOOD  01 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
-O.OOOOOOOOD  00 
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XB<  5)  *  0.95514D  00 

XB(  6)  «  0.25776D  01 

XB<  7)  *  0.60578D  00 

XB<  8)  «  0.32747D  01 

X8<  9)  *  0.33796D  00 

THE  CURRENT  VALUE  Of  THE  OBJECTIVE  FUNCTION  *  -0.62621B  01 

THE  REDUCED  COSTS 

THE  1-TH  REDUCED  COST  *  -0.28866D-14 

THE  2-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  3-TH  REDUCED  COST  *  -O.OOOOOD  00 

THE  4-TH  REDUCED  COST  =  0.28866D-14 

THE  5-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  &-TH  REDUCED  COST  =  -0.20000D  01 

THE  7-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  8-T  , .EDUCED  COST  =  -O.OOOOOD  00 

THE  9-iH  REDUCED  COST  =  -O.OOOOOD  00 

THE  10-TH  REDUCED  COST  =  -0.20000D  01 

THE  11-TH  REDUCED  COST  -O.OOOOOD  00 
THE  12-TH  REDUCED  COST  =  -0.20000D  01 

THE  13-TH  REDUCED  COST  =  -0.20000D  01 

THE  14-TH  REDUCED  COST  =  -0.20000D  01 

THE  15- TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  14-TH  REDUCED  COST  =  -O.OOOOOD  00 

THE  CURRENT  BASIC  SOLUTION  IS  FEASIBLE  AND  HENCE  OFTIHAL. 

THE  NONZERO  VARIABLES  ARE  AS  FOLLOWS: 


X( 

0)  = 

-0.62621D  01 

X(  50)  = 

-O.OOOOOD  00 

X( 

5)  = 

0.20196D  00 

X<  11)  = 

0.33537D-01 

X< 

7)  * 

0.17223D  00 

X( 

3)  * 

0.95514D  00 

X( 

9)  * 

0.25774D  01 

X( 

2)  = 

0.60578D  00 

X( 

8)  = 

0.32767D  01 

x< 

15)  * 

0.33796D  00 

THE  OPTIMAL  VALUE  OF  THE  OBJECTIVE  FUNCTION  IS  -0.62621D  01 
AUXILIARY  RESULTS  OF  THE  MINI  MONTE  CARLO  STUDY 


VALUES  OF  DELTA  BETA  STAR 


SAMPLE  NUMBER 

= 

1 

-0.01197 

-0.20926 

SAMPLE  NUMBER 

3 

2 

0.51131 

-0.10930 

SAMPLE  NUMBER 

5 

3 

1.43557 

-0.71028 

SAMPLE  NUMBER 

3 

4 

0.61193 

-0.44646 

SAMPLE  NUMBER 

3 

5 

0.61553 

-0.51802 

SAMPLE  NUMBER 

3 

6 

-0.28073 

-0.21389 

SAMPLE  NUMBER 

3 

7 

-0.18772 

-0.34817 

SAMPLE  NUMBER 

S 

8 

-0.36551 

-0.05786 

SAMPLE  NUMBER 

3 

9 

0.95646 

-0.28508 

SAMPLE  NUMBER 

3 

10 

0.03818 

-0,28217 

SAMPLE  NUMBER 

3 

11 

1,59676 

-0.62135 

SAMPLE  NUMBER 

3 

12 

0.62508 

-0.34275 

SAMPLE  NUMBER 

3 

13 

-0.10425 

-0.41664 

SAMPLE  NUMBER 

3 

14 

0.67201 

-0.34557 

SAMPLE  NUMBER 

3 

15 

1.18308 

-0.61111 

SAMPLE  NUMBER 

3 

16 

1.24476 

-0,56381 

SAMPLE  NUMBER 

3 

17 

0.66750 

-0,66472 

SAMPLE  NUMBER 

3 

18 

1.85856 

-0.79734 

SAMPLE  NUMBER 

3 

19 

0.99733 

-0,72490 

SAMPLE  NUMBER 

s 

20 

1.12891 

-0.81330 

ESTIMATED  COVARIANCE  OF  DELTA  BETA  STAR 
0.448991  -0.122225 

-0.122225  0.054857 

SUM  OF  THE  OPTIMAL  OBJECTIVE  FUNCTIONS  OVER  ALL  SAMPLES  = 
MAIN  RESULTS  OF  THE  MINI  MONTE  CARLO  STUDY 

ESTIMATED  VALUE  OF  SIGMA  *  2.817978 

AUXILIARY  RESULT  5  THE  ESTIMATE  OF  SIGMA  5  *  7.736290 

ESTIMATED  COVARIANCE  OF  THE  REGRESSION  PARAMETER  VECTOR.  BETA 
3.565439  -0.970588 

-0.970588  0.435623 


131,427570 
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